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ABSTRACT 

We present an advance in the u se of C assini observations of stellar occultations by the rings 
of Saturn for stellar studies. Stewart et al. ( |2013| ) demonstrated the potential use of such ob¬ 
servations for measuring stellar angular diameters. Here, we use these same observations, and 
tomographic imaging reconstruction techniques, to produce two dimensional images of com¬ 
plex stellar systems. We detail the determination of the basic observational reference frame. A 
technique for recovering model-independent brightness profiles for data from each occulting 
edge is discussed, along with the tomographic combination of these profiles to build an im¬ 
age of the source star. Finally we demonstrate the technique with recovered images of the a 
Centauri binary system and the circumstellar environment of the evolved late-type giant star, 
Mira. 
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1 INTRODUCTION 

Often, a great deal of useful information about celestial targets may 
be obtained from sparse imaging data, such as one-dimensional 
brightness profiles. A novel example of such data was presented 
in [Stewart et al.| ( |2013f (Paper I) which demonstrated the use of 
Cassini-Y IMS (Visible and Infrared Mapping Spectrometer) obser¬ 
vations of stellar occultations by the Saturnian rings to recover high 
spatial resolution information about the circumstellar environment 
of Mira. 

However, there are circumstances in commonly encountered 
astronomical systems in which simple assumptions such as spher¬ 
ical symmetry of a target object are no longer reliable. In some 
of the most extreme examples, such as imaging the close environ¬ 
ments of very dusty mass-losing evolved stars ( [Tuthill et al.|2005| ), 
the recovered morphology may have very few residual symmetries 
or predictable structural elements with which to build models for 
low-order parameter fitting. In such cases, the only recourse is some 
form of model-independent image recovery process. 

Although Paper 1 primarily employed the Cassini occultation 
data to recover a simple one-dimensional brightness profile, the 
data themselves do explore more of the object’s flux distribution 
than this. In particular, over the course of single passage of the en¬ 
tire Saturnian ring system across the sightline to a celestial target, 
we found that there were usually a handful of usable occultation 
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events at various known sharp edges within the ring structure (see 
Paper 1 for an explanation of the basic occultation geometry). Most 
importantly, the apparent sky-orientation of the occulting edge for 
each occultation event is a different angle, so that an ensemble of 
such events - built up over the course of one or several complete 
passes behind the rings - yields a data set capable of constraining 
a full two dimensional image. 

More generally the process of reconstructing an image from a 
series of overlapping one-dimensional projections is known as To¬ 
mography. This has been particularly important in medical imaging 
for the past several decades as a non invasive way to view structures 
internal to the human body. It turns out that the data recovered from 
Cassini-Y IMS occultations form a sufficiently similar representa¬ 
tion that we are able to employ these same techniques to derive 
two dimensional images of stellar sources and their circumstellar 
environments. 

The first problem confronted in the image formation process 
is recovery of the angle on-sky from each occulting edge in a suit¬ 
able celestial reference frame. This is discussed in Section [2] The 
process of generating tomographic reconstructions of stellar targets 
themselves is detailed in Section [3] and is broken into steps. The 
first (Section[3T} is to recover the one-dimensional brightness pro¬ 
file of the source star from each occultation observation. These are 
then combined tomographically to build an image (Section [T2] ) in 
an iterative process to determine the projection registration (Sec¬ 
tion [3T3j. 
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Figure 1. An example of the two coordinate systems used to compare vec¬ 
tors from different sources. The blue vectors identify the f,g,h coordinate 
system, with h along the line of sight from Cassini to the star, and centred 
on the point where this line intercepts the ring plane, the RPI (Ring Plane 
Intercept). The f,g coordinates define Cassini’s sky-plane with / lying in 
the ring plane. The black vectors show the x,y,z coordinate system with its 
origin at the centre of the planet, z through the poles and x defined by the 
ring plane projection of the vector towards the star. The red line (r) joins 
the origins of the two coordinate systems and 0 is the angle between the 
x axis and the kronocentric longitude of the RPI. The red dashed vector is 
the standard kronocentric x direction in the ring plane, and 0* is the angle 
between this and the kronocentric longitude of the star. The projection of 
the h axis into the ring plane is parallel to the x axis, and the angle between 
the h axis and this projection is B *. 


2 POSITION ANGLE OF PROJECTIONS 



Figure 2. Identifying the important vectors in Cassini’s sky-plane (defined 
by the / and g axes), g is the projection of Saturn’s north vector translated to 
the RPI, where the line of sight to the star (unlabelled h axis) is intercepted 
by the planetary ring-plane. is the direction towards the centre of Saturn. 
v_l is the direction normal to the rings which corresponds to the direction of 
spatial information encoded in the occultation N is the terrestrial celestial 
north vector translated to the RPI and P.A. is the position angle of v± relative 
to the Earth’s north. Base image: PIA12518 ( jCassini Imaging Team|2010| 

onal axes (/ and g) in Cassini's sky-plane. This coordinate system 
and the relevant vectors are shown in Figure [2] 

The sky-plane coordinates can be defined in terms of the Kro¬ 
nocentric coordinates by 


In order to make sense of measured asymmetries in fitted bright¬ 
ness profiles, and to register structure against published literature, 
the projection angle of each one-dimensional profile must be deter¬ 
mined. This is the angle normal to the occulting edge as viewed in 
the sky-plane of Cassini , relative to north in the celestial sphere 
of the Earth. The spacecraft position and pointing direction are 
recorded by mission controllers in a kronocentric reference frame, 
so the terrestrial celestial north vector must be determined relative 
to that of Saturn’s north. 

We start by defining two coordinate systems: a kronocentric 
reference frame with its origin at the centre of Saturn, and another 
centred on the point where the planetary ring plane intercepts the 
line of sight from Cassini to the star as shown in Figure [T] Using 
these two coordinate systems simplifies the process of determining 
the required angle using vectors coming from different sources. 

The kronocentric coordinate system has its origin at Saturn’s 
centre (f} c ) with z through the planet’s poles and x and y in the 
ring plane. The standard 0° longitude in the kronocentric coordi¬ 
nate frame is defined to be the ascending node of the ring-plane on 
the terrestrial equatorial plane in J2000. In order to simplify the ge¬ 
ometry, we instead define the x axis to be the ring-plane projection 
of the vector from toward the star. This direction is effectively 
fixed for any particular star, varying only with planetary preces¬ 
sion and stellar proper motion on million year time-scales. This 
is achieved by rotating by the kronian longitude of the star (0*) 
around the z axis, 


f = y-y r (2) 

g = zcos(2?*) — (x — x r ) sin(Z?*) (3) 

where x r and y r are the x, y coordinates of the RPI and B * 
is the kronocentric latitude of the star which is constant for any 
particular star. This is effectively just a rotation of — B * around the 
y axis as 

[hJ,g]=^ 2 {-B*)[x-x r ,y-y r ,z] (4) 

The determination of the radius (r) and longitude (0) of the 
of the RPI is discussed in |Stewart et al.| ( |2013| ). From these, the 
determination of x r and y r is trivial. 

r= (rcos(0),rsin(0),O) = (x,y,z) r (5) 

We are now able to rotate the terrestrial celestial north vector 
into Cassini's sky-plane by employing three matrix rotations. Two 
of the three matrices required are taken from Equation [l] and Equa¬ 
tion [4] while the final matrix is C eP((X p , 8 p ) which rotates J2000 
coordinates into kronocentric equatorial coordinates. These trans¬ 
formations may be combined into a single matrix operation: 

( 6 ) 

enabling translation of the Earth’s celestial north into the 
RPI/sky-plane coordinate system. 

The position angle (P.A.) of the occultation can be determined 
as follows, 


[x,y,z] = K. 3 (0*)[*,;y,zJ G) 

The second coordinate system has its origin at the ring-plane 
intercept (RPI) which is the point where the planetary ring plane is 
intercepted by the line of sight from Cassini to the star, with one 
axis (h) running along the line of sight to the star, and the orthog¬ 


P.A. = 


cos 1 (v_|_ • N) 
v ± 


(7) 


where is defined to be the direction normal to the projection 
of the ring edge into the sky plane and is different for each occulted 
edge. This can be found by 
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Figure 3. A comparison of angular diversity for representative occultation 
geometries using indicative edges. The long red vectors indicate the appar¬ 
ent trajectory of the star behind the rings. The short cyan vectors identify 
the normal to some of the occulting edges, (a) shows a minimal occultation 
with only two edges that have almost parallel PA. s. (b) shows an almost 
identical orbital geometry, with vastly different angular diversity, although 
most edges produce almost redundant P.A. s. The cyan vector clusters to 
the top left of each shows an angular diversity rose, a quick visual method 
introduced here to compare angular diversity between observations. Base 
image: PI A10446 ( [Cassini Imaging Team|2010) 


v ± = [0,sin(vO,cos(vO] = l h J,g]v± ( 8 ) 

where yt is determined from 

tan (V) = -^I (9) 

af 

The movement of the line of sight to the star in the ( f,g ) plane 
is 

/ = — Arsin(0) + A0rcos(0) (10) 

g= [Arcos(0) + A0rsin(0)] sin(2?*) (11) 

with Ar and A0 being the change in position of the RPI. Ar 
is locally constant whilst A0 varies across the ring edge. Differen¬ 
tiating Equations [T0| and |TT| with respect to r and combining into 
Equation [9]produce s 

d,g 

tan(yr) = -— = -tan(0)sin(5*) (12) 

cij 

which can be used to solve Equation[8]to produce for each 
occulting edge. 

This set of transformations results in the h axis lying in direc¬ 
tion of the line of sight to the star, and so it may be ignored here. 
With and N now rotated to lie in the sky-plane ( f,g ), obtaining 
the angle between them is straightforward. The convention adopted 
is that the angle will be 0° if they happen to align; it will be pos¬ 
itive if is to the left of N toward the east in the sky, otherwise 
negative (toward the west). 


2.1 Angular Diversity of Observations 

As the star is observed to pass behind the rings, it encounters many 
edges. Each of these edges yields an occultation lightcurve from 
which a brightness profile can be obtained with a different P.A., 
thus sampling the star in different directions. Every transit behind 
the ring plane will produce a different number and diversity of an¬ 
gles. Figure [3] shows an example of the differences in angular di¬ 
versity between two occultations with similar orbital geometry. The 
occultation on the left, labelled (a), crosses only two sharp edges, 
and produces an angular diversity rose with poor coverage, and al¬ 
most no sampling in the N-S direction. The vector direction in the 
“rose” is the direction of spatial sampling of the stellar image which 
is orthogonal to the projection direction. Occultation (b) has much 
greater angular diversity, although many of the angles are largely 
redundant. It is has far greater sampling in the N-S direction than 
it has E-W, and is generally superior to occultation (a) in almost all 
respects for image recovery. 


3 TOMOGRAPHIC IMAGE RECONSTRUCTION 

3.1 Model-Independent Brightness Profile Recovery 

Occultations typically recover one-dimensional brightness profiles. 
These are the projection of the two-dimensional image of the star in 
the sky-plane onto a particular axis. Paper I demonstrated the suc¬ 
cessful recovery of circularly symmetric model brightness profiles 
from Cassini occultation observations. The use of models is a par¬ 
ticularly potent way to exploit data which are sparse, noisy or both. 
Best-fit models are able to contribute significantly to studies of tar¬ 
get structure, but they are only capable of parameterising structural 
elements that are already foreseen in the model construction. 

On the other hand, the present study attempts to obtain a 
“true” tomographic image, without enforcing symmetries and low- 
order representation that invariably come with the fitting of mod¬ 
els. As a first step, we need to recover one-dimensional brightness 
profiles for each occulting edge, independent of any model. The 
ideal way to perform this would be to deconvolve the observed 
lightcurve with the point source lightcurve. Unfortunately the de- 
convolution process produces non-unique solutions and is severely 
under-constrained due to the presence of even small amounts of 
noise found in all measurements. For this reason, instead of per¬ 
forming a relatively straightforward Fourier deconvolution, more 
sophisticated algorithms are required to recover model-independent 
brightness profiles. 

The approach taken is to allow a model brightness profile to 
include a free parameter for each data point in the observation. 
Each point is allowed to vary between zero and one before be¬ 
ing convolved with the point source lightcurve to produce a model 
lightcurve. The point source lightcurve is determined using Fresnel 
diffraction and the known orbital geometry of Cassini as detailed 
in [Stewart et al.||2013) . These model lightcurves are in turn com¬ 
pared against the observed lightcurve. The process iterates by al¬ 
lowing the synthetic brightness profile to evolve towards a better 
fit. This ultimately produces a brightness profile which generates 
a lightcurve closely matching that which has been observed. Some 
examples of this process are shown in Figure]?] 

The model fitting was performed with the minimize function in 
Scipy’s optimize module. This is a series of maximum likelihood fit¬ 
ting algorithms which are a native part of the Scipy scientific com¬ 
puting environment for the Python language. The ability to regener¬ 
ate accurate brightness profiles using this method was exhaustively 
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Simulated Point Source Recovered 

Lightcurves Lightcurve Brightness Profiles 



Figure 4. Examples of the model-independent recovery of brightness pro¬ 
files (right), from simulated lightcurves (left) and the corresponding point 
source lightcurve (centre). The recovered profiles are used in Figure [5] and 
the simulated lightcurves are produced from the ’Original Image’ in that 
figure, in the directions indicated by the angular diversity rose. The position 
angle of each lightcurve is labelled on the left hand side in the correspond¬ 
ing colour. The occultation geometry used for each lightcurve simulation is 
identical, resulting in a common point source lightcurve. 

tested by generating synthetic brightness profiles, adding noise, and 
fitting. At a signal-to-noise ratio of 100 we are able to recover the 
overall structure of the brightness profile. When the signal-to-noise 
ratio is reduced to 10, we are able to recover the envelope of the 
brightness profile, but not finer details. 

3.1.1 The Point Source Case 

Ideally an experimental test of the instrumental transfer function 
would be performed in order to validate the model used to trans¬ 
late observational data into measurements of target morphology; in 
particular to confirm that departures from idealised representations 
of the spacecraft detection system and ring edges do not adversely 
affect the data. Often, observations of a known unresolved point 
source reference star are taken to serve this purpose. However, the 
only stars bright enough to be observed with VIMS in occultation 
mode, are inherently large enough to be resolved, making a true 
point source observation impossible. This is a problem sometimes 
faced by long-baseline ground-based interferometers, where it can 
be difficult to identify a calibrator with sufficient surface bright¬ 
ness to be detected whilst remaining unresolved. In the absence of 
a true point-source, the most rigorous test available is a comparison 
of fitted model parameters to literature values over a wide spec¬ 
tral range. This was performed by comparing uniform disc fits to 
terrestrial interferometric observations of Mira and was presented 
in [Stewart et al.| ( [2013| ). Close agreement with prior literature was 
shown to validate results produced from occultations. 


3.2 Back Projection and Regression 

The building of a tomographic image requires two steps. These are 
back projection and regression. Back projection is projecting the 
recovered brightness profiles over each other in the direction de¬ 
termined by the P.A. The overlapping profiles immediately form 
a crude image of the object. More projections and increased an¬ 
gular diversity provides a visibly better image at this stage, but a 
limited number of projections or poor angular diversity, as always 
occurs with these observations, makes the regression step crucial. 
Regression allows the image to be modified in some way, such as 
smoothing or concentrating, whilst maintaining the fit to the back 
projections. 

Tomography is a mature field and has many existing algo¬ 
rithms and tools capable of recovering images from well sampled 
data. Most of these tools are not suited to the sparse angular di¬ 
versity these Cassini observations produce. The selection and im¬ 
plementation of a suitable regression algorithm is therefore criti¬ 
cal. The Scikit-learn project ( [Pedregosa et al.|20TT] > offers machine 
learning regression algorithms in Python which proved to be suf¬ 
ficiently versatile and capable for image reconstruction. Their im¬ 
plementation of an elastic net linear regularised regression algo¬ 
rithm ( [Hui Zou||2005) ) was found to be highly effective with the 
sparse data produced with this technique. The selection of this reg- 
ulariser was made after extensive testing of reconstructions of syn¬ 
thetic objects and found visually to provide the fewest artefacts and 
the best defined real structure. 


3.3 Centroiding Brightness Profiles 

The brightness profiles as recovered by the algorithms discussed 
in Section [3TT| possess no unambiguous phase reference, with the 
exception of simple stellar geometries such as binary systems, 
or simple circularly symmetric stellar targets. That is to say that 
when attempting to back-project more than two brightness profiles 
taken at different angles to synthesise a two-dimensional image, 
the phase registration or shift of each projection is not fixed, a sce¬ 
nario not commonly encountered in the tomography literature. For 
more complex stellar environments the image reconstruction re¬ 
quires some sort of referencing between projections, otherwise se¬ 
vere artefacts can be introduced in the resultant image. If the object 
being imaged is known to be inversion symmetric, this process is 
quite straightforward and can be accomplished in several different 
ways with relative ease. These include centroiding based on fitting 
Gaussian functions, or registration based on the centre of mass of 
each of the profiles. 

For complex objects known to be dominated by significant 
asymmetries, or those with unknown geometries, this process can 
result in misfits and and introduction of artefacts. The exact regis¬ 
tration of an ensemble of angular-diverse brightness profiles cannot 
be unambiguously determined for an unknown object. In order to 
assess the impact of registration on the quality of the output images, 
a study simulating the entire process from occultation lightcurve to 
image was performed. 

Results from these simulations are given in Figure [5] which 
depicts the impact upon the final image reconstruction of differ¬ 
ing levels of centroiding accuracy. In this case, a small picture of 
Saturn (top) was adopted as a useful illustrative model image as 
it exhibits a distinct asymmetric geometry, with both coarse and 
fine features ideally suited to benchmarking image recovery tech¬ 
niques. The “True” centroiding aligns the profiles based on exact 
foreknowledge of the image (and is therefore only possible in a 















Stellar Imaging with Cassini Occultations II 5 



Sinograms 


-50 



'True' Centroiding 


120 mas 




Modest Misalignment 



120 mas 


Iterative Centroiding 


120 mas 




- -50 

0 40 80 120 

Projection Angle (°) 



Figure 5. Simulated images demonstrating the impact of centroiding on 
image recovery. The top frame shows the original model image that was 
used to generate one-dimensional projections. The angular diversity rose 
in the top right of that frame shows the five projection directions. These 
projections were convolved with a point source light curve produced from 
a typical occultation geometry as shown in Figure [4] Realistic noise was 
added, which remains visible in the reconstructed images as high frequency 
fluctuations. These simulated lightcurves were then fitted with the model- 
independent process detailed in Section [TT| to produce brightness profiles. 
The brightness profiles recovered for each projection angle can be seen in 
the sinograms down the right hand side with corresponding recovered im¬ 
ages to their left. The vertical axis on the sinograms is at the same scale as 
the recovered images. 


simulation), and shows the upper bound of expected image quality 
with perfect registration. The Gaussian centroiding process fits a 
Gaussian function to the recovered brightness profile and performs 
a registration based on the function centre. For this particular case 
(and any significantly asymmetric object), this approach is not op¬ 
timal as the brightest part of the image does not occur at the cen¬ 


tre and as a consequence the image exhibits artefacts. For example 
the ring structure recovered is now asymmetrically disposed with 
respect to the planet’s disc. The third image/sinogram pair demon¬ 
strates that even modest misalignment of the registration can pro¬ 
duce artefacts (for example, loss of detail like the elliptical shape 
of the rings). 

After trying several strategies, the most effective method for 
registering test images was determined to be an iterative process 
using recovered tomographic images as discussed in Section |3l2| 
At the first iteration, Gaussian centroiding (as above) delivered a 
trial recovered image. Next a filtered image was produced by ap¬ 
plication of a two-pixel Gaussian blur (removing fine detail), and 
a binary mask applied at the 15% cut level truncated low level pe¬ 
ripheral noise and artefacts. At the fitting step, the registration of 
each recovered one-dimensional profile was treated as a free pa¬ 
rameter used to minimise the misfit with the filtered image. Based 
on this new ensemble of registrations, a new tomographic image 
was recovered and the algorithm had progressed a full cycle. This 
whole process was iterated until none of the profiles were able to 
provide a better fit with a change in offset. Given the existence of a 
known ‘truth’ image in our simulations, we were able to show that 
the algorithm converged to near the known correct registrations in 
only a few iterations. The results of this process are shown in the 
bottom row of Figure [5] where the overall shape and balance of the 
image is recovered with minimal artefacts. 


4 RESULTS OF KRONOCYCLIC TOMOGRAPHY 
4.1 Binary Systems: a Centauri 

The a Centauri system is the closest stellar system to our solar 
system. It appears to the eye as the third brightest star in the sky, 
and the binary nature of a Centauri AB has been observed for over 
three hundred years. The orbital elements of the system have been 
well defined as it is easily resolved with a small telescope and has 
a long history of measurements. Modern techniques have been able 
to substantially refine this, providing highly accurate position pre¬ 
dictions. 

The effectiveness of tomographic imaging from Cassini stel¬ 
lar occultations was tested by comparing reconstructions from ob¬ 
servations of a Centauri at two epochs to the orbital predictions of 
|Pourbaix et al.| ( |2002] ). The first of these reconstructions is shown in 
Figure[6] from an occultation observed on 29 April 2008 ( Cassini’s 
revolution 066). The second is from 6 March 2009 ( Cassini's revo¬ 
lution 105) and is shown in Figure [7] These reconstructions found 
the pair of stars to be located as expected to within 5% in separa¬ 
tion, and 2° in position angle as shown in Table [I] 

The a Centauri image reconstructions are particularly hin¬ 
dered by the limited angular diversity provided by the specific oc¬ 
cultation geometry. Both epochs suffered from a lack of diversity 
amongst the projection angles, limiting the ability of the technique 
to accurately reconstruct the correct position angle. For R. 066 the 
projections were all within a 25° angle, and for R. 105 they were 
within 8.1° as shown by the angular diversity roses in Figures [6] 
an d0 Both cases also had all their projection angles within 19° of 
the binary position angle, increasing the uncertainty on the separa¬ 
tion. This would be improved if the position angle and projection 
angles were closer to perpendicular. We expect that improved angu¬ 
lar diversity would provide substantially closer agreement with the 
predicted binary separation and position angle. In spite of these lim¬ 
itations, our observations at these two epochs reveal the expected 
orbital progression of the binary system. 
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Figure 6. A tomographic reconstruction of the a Cen system from obser¬ 
vations performed on 29 April 2008 during Cassini’s revolution 066. The 
axes are centred on a Cen A, and expressed in arcseconds of right ascension 
and declination, with north up, and east to the left. The red pluses indicate 
the relative position of a Cen B at the start of each calendar year based on 
an ephemeris from |Pourbaix et al .1(2002) . The blue cross, labelled with the 
decimal year, shows the predicted relative position of a Cen B at this epoch. 
The insets show a 1" zoomed view of each of the stars and the tomographic 
artefacts in their vicinity. 



Figure 7. A tomographic reconstruction of the a Cen system from data 
taken on 6 March 2009 during Cassini’s revolution 105. All details as per 


Figure[6] 





Epoch 

Separation (") 

P.A. ( 

°) 

(Year) 

Predicted 

Measured 

Predicted 

Measured 

2008.37 

8.007±0.007 

7.71±0.83 

238.55±0.09 

236.5±2.9 

2009.18 

7.399±0.007 

7.79±0.55 

241.60±0.11 

243.6±2.4 


Table 1. Predictions and results of binary parameters from Alpha Centauri 
comparing the |Pourbaix et al. 112002) ephemeris and tomographic images. 


4.2 Evolved Stars: o Ceti 


o Ceti is the archetype of the Mira Variable stellar class. These 
are evolved stars which have exhausted both the hydrogen and he¬ 
lium in their cores, and are known to exhibit complex circumstellar 
structure ( |Tuthill et al,12005| ). Mira itself has been shown to have 
an asymmetric circumstellar environment (Karovska et al.||1997[ 
Martin et al. |2007|) inc luding dusty concentric shell-like struc¬ 
ture s ([Lopez et al.|1997[|Cruzalebes et al.|1998| . 


|Stewart et al.| ( [2013] ) demonstrated the use of Cassini stellar 
occultations to measure the wavelength dependant mean diameter 
of Mira. Here we use the same observations from Cassini's revolu¬ 
tion 10 to tomographically reconstruct a high resolution image of 
Mira. Observations from earlier epochs (Revs. 8 and 9) were found 
to each have relatively poor angular diversity. The inherent vari¬ 
able nature of Mira prevents using brightness profiles from differ¬ 
ent epochs being used together to improve this diversity. The results 
of this reconstruction are shown in Figure [8] revealing the complex 
structure of the circumstellar environment of Mira. Occultations by 
the outer edge of the Keeler Gap on ingress, and the Encke Gap on 
both ingress and egress were used for this reconstruction. These 
particular edges were found to produce a consistent reconstruc¬ 
tion with minimal noise. It shows a bright core with a diameter 
of around 30 mas as shown in |Stewart et al.| ( [T013| > which contain 
the vast majority of the stellar flux. The remaining flux forms a 
complex dusty environment cooling in discrete steps as a function 
of radius. The first of these occurs at a radius of around 80 mas and 
the second at a radius of around 200 mas. These are consistent with 
previously identified dynamic H 2 O shell structures ( [Matsuura et al.| 
|2002| ). All of these structures are revealed to be somewhat asym¬ 
metric with their respective radii varying with angle. The most ob¬ 
vious of these asymmetries is the enlargement of the outer cooler 
layer to the north-east of the star. Similar asymmetries within the 
inner 100 mas have previously been identified in the visible part of 
the spectrum ( [Wilson et al.|1992[|Karovska et al.|1997) . Figure [8] 
also shows examples of some of the artefacts which this process 
can introduce. One of these is the dumpiness of the dusty areas on 
approximately the same scale as the angular sampling. Another is 
the rays which can be seen stretching across the image perpendicu¬ 
lar to the sampling directions indicated in the angular diversity rose. 
These are a by-product of back projection, and are reduced with im¬ 
proved angular diversity. The consequence of these artefacts is that 
fine structure in the resulting images should be treated with scepti¬ 
cism, however the overall envelope appears to be robust. A detailed 
astrophysical interpretation, including polychromatic reconstructed 
images, will be presented in a forthcoming paper. 


5 CONCLUSIONS 

We have successfully demonstrated a technique to use stellar occul¬ 
tations by the rings of Saturn, as observed by Cassini , to produce 
tomographic images of stellar systems. This required a thorough 
understanding of the spacecraft’s orbital geometry, and a transla¬ 
tion into a terrestrial reference frame. It also entailed the extraction 
of model-independent one dimensional brightness profiles and de¬ 
veloping an understanding of tomographic processes. The process 
was tested and verified by imaging the known geometry of the a 
Centauri binary system. It was then applied to the dynamic, ever 
evolving inner regions of Mira’s dusty circumstellar environment. 

This technique has the potential to obtain polychromatic stel¬ 
lar images from a space platform at similar angular resolutions to 
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Figure 8. A three colour tomographic reconstruction of Mira from 29 June 2005 during Cassini’s revolution 10. Red is 4.66 fim, green is 3.32 fim, and blue 
is 1.34 j am. North is up and east is left and the axes are in milliarcseconds of right ascension and declination. The angular diversity rose in the top right of the 
image shows the projection directions used in the reconstruction. The sampling resolution of each projection is indicated by the length of the corresponding 
stroke in the rose. 


terrestrial long-baseline interferometry. It can be applied to both 
archived and new observations (until Cassini’s planned demise in 
2017) providing they have appropriate occultation geometry. Its ap¬ 
plication to other stellar systems could potentially yield reveal ad¬ 
ditional spatial information about circumstellar environments. 
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